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Abstract 

We show that Pyragas delayed feedback control can stabilize an 
unstable periodic orbit (UPO) that arises from a generic subcritical 
Hopf bifurcation of a stable equilibrium in an n-dimensional dynam- 
ical system. This extends results of Fiedler et al. [PRL 98, 114101 
(2007)], who demonstrated that such feedback control can stabilize 
the UPO associated with a two-dimensional subcritical Hopf normal 
form. Pyragas feedback requires an appropriate choice of a feedback 
gain matrix for stabilization, as well as knowledge of the period of the 
targeted UPO. We apply feedback in the directions tangent to the two- 
dimensional center manifold. We parameterize the feedback gain by 
a modulus and a phase angle, and give explicit formulae for choosing 
these two parameters given the period of the UPO in a neighborhood 
of the bifurcation point. We show, first heuristically, and then rigor- 
ously by a center manifold reduction for delay differential equations, 
that the stabilization mechanism involves a highly degenerate Hopf bi- 
furcation problem that is induced by the time-delayed feedback. When 
the feedback gain modulus reaches a threshold for stabilization, both 
of the genericity assumptions associated with a two-dimensional Hopf 



bifurcation are violated: the eigenvalues of the linearized problem do 
not cross the imaginary axis as the bifurcation parameter is varied, 
and the real part of the cubic coefficient of the normal form vanishes. 
Our analysis of this degenerate bifurcation problem reveals two qual- 
itatively distinct cases when unfolded in a two-parameter plane. In 
each case, Pyragas-type feedback successfully stabilizes the branch of 
small-amplitude UPOs in a neighborhood of the original bifurcation 
point, provided that the phase angle satisfies a certain restriction. 



1 Introduction 

Chaotic attractors typically possess a dense set of unstable periodic orbits 
(UPOs). This form of phase space skeleton was exploited in a control scheme 
developed by Ott, Grebogi and Yorke [I] in the 1990s. Their approach pro- 
vided a method for stabilizing targeted UPOs of chaotic systems for ap- 
plication in both numerical simulations and laboratory experiments. This 
method spawned the development of a number of related and alternative con- 
trol schemes, with similar goals, which can be utilized in systems for which 
strictly periodic behavior is attractive |2]. One of these schemes, which has 
been especially well investigated and tested, was proposed by Pyragas [3]. 

Pyragas control exploits the symmetry of a periodic orbit in a natural 
way by providing, in its simplest realization, additive feedback in the form 



K 



x(t - r) - x(t) 



Here x(t) G M" is the state vector of the dynamical system at time t, t is 
the period of the targeted UPO, and K is a constant n x n feedback gain 
matrix. The scheme is manifestly noninvasive, since the feedback vanishes 
when the system reaches the r-periodic target state. Setting aside the dif- 
ficult questions related to basins of attraction, there are then just two key 
ingredients to the successful implementation of this approach: the period r 
of the targeted UPO is needed, and the feedback gain matrix K needs to 
guarantee stabilization. Only to the extent that an appropriate choice of K 
is required does the method rely on detailed knowledge, beyond the period, of 
the structure of the UPO in phase space. For a review of the extensive liter- 
ature on applications of Pyragas feedback, including successful experimental 
implementations, see 
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This paper is motivated by the question of how to choose the feedback gain 
K in Pyragas control to ensure that it will be effective. We focus on a simple, 
generic mechanism for the creation of an unstable periodic orbit: the sub- 
critical Hopf bifurcation of a stable equilibrium. Other generic mechanisms 
for creating UPOs in dynamical systems include homoclinic bifurcations, 
saddle-node (or fold) bifurcations of limit cycles, saddle-node bifurcations 
of fixed points on an invariant circle, and period-doubling bifurcations [5]. 
There have been a number of successful demonstrations of Pyragas control 
of periodic orbits destabilized through a period-doubhng bifurcation (see, for 
instance, [6]). Postlethwaite has shown that Pyragas- type feedback can 
stabilize a UPO arising from a subcritical bifurcation from a robust hetero- 
clinic cycle in a three-dimensional system of equivariant ordinary differential 
equations. Interestingly, Pyragas feedback works in this case even though 
the period of the targeted orbit, and hence the time-delay, diverges as the 
heteroclinic bifurcation point is approached. Fiedler et al. [8j have inves- 
tigated an example of Pyragas control that stabilizes a circular limit cycle 
{i.e. a "rotating wave") near a fold bifurcation in a planar system of ordi- 
nary differential equations with 0(2)-symmetry, and successfully applied this 
to a higher-dimensional model taken from nonlinear optics that possesses a 
similar rotational symmetry. 

The first example that demonstrated the successful stabilization of a UPO 
arising from a subcritical Hopf bifurcation was given in [S], and further an- 
alyzed in [TU], with an experimental implementation described in [TT]. In 
these papers, the authors added Pyragas feedback directly to the Hopf nor- 
mal form: 

^ = (A + i)z{t) + (1 + ii)\z{t)\^z{t) + 6oe^^ [{z{t - r) - z{t)\ . (1) 

Here r = ^^^^ is the period of the UPO, and the complex number boe^^^ 
plays the role of the feedback gain matrix K. A beauty of this simple ex- 
ample is that it represents a rare instance in which solutions of a nonlinear 
delay differential equation can be computed analytically in closed form, and 
their bifurcations can be studied with comparable finesse. Specifically, us- 
ing methods of bifurcation theory, the authors were able to understand the 
mechanism for stabilization in this example. For instance, they showed that 
the feedback control leads to additional delay-induced Hopf bifurcations of 
the equilibrium z = 0, and consequently it is possible to change the equilib- 
rium's stability so that the original subcritical Hopf bifurcation to the UPO 
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turns into a supercritical bifurcation to a stable periodic orbit. An impor- 
tant contribution of the Fiedler et al. [9J paper was that it also provided a 
counterexample to a published claim \l2l that Pyragas control is impossible 
when the UPO has an odd number of real positive Floquet multipliers greater 
than one. In the ten years between the published claims of the odd number 
limitation [12] and the first counterexample to it [9], a number of modifica- 
tions of the Pyragas control scheme were developed. One of these, based on 
introducing an additional unstable direction via the controller, was proposed 
in order to stabilize UPOs created by a subcritical Hopf bifurcation [131 E]) 
including an example applied to the Lorenz equations [T5] . 

Recently an analysis of the Lorenz equations with the standard Pyragas 
feedback provided a second example of stabilization of a UPO resulting from 
a subcritical Hopf bifurcation. Specifically, Postlethwaite and Silber [16] 
demonstrated that the stabilization mechanism identified by Fiedler et al. [9] 
can also apply to UPOs in higher- dimensional systems, provided that the 
feedback gain matrix is chosen correctly. The strategy they outlined is to 
add feedback of the type investigated in [HI HO] in the directions tangent to 
the center manifold of the uncontrolled Lorenz system. Stabilization of the 
UPOs is then possible over a broad range of control parameter values. 

The reduction of higher- dimensional systems to the two-dimensional nor- 
mal form near a Hopf bifurcation is a standard procedure [17] . Likewise, such 
systems with additive Pyragas feedback, now infinite-dimensional, can also 
be reduced to the standard two-dimensional normal form in the vicinity of a 
Hopf bifurcation, where the parameters of the feedback gain matrix modify 
the coefficients in the normal form. In the Fiedler example [9], the feedback 
terms are added directly to the Hopf normal form, but a surprising result 
of [in] was that the same sequence of bifurcations identified in the simpler 
normal form example also appears in this higher-dimensional example. This 
result is generalized further in the current paper. 

A further motivation for the work we present in this paper is to under- 
stand the origin of a particular degenerate Hopf bifurcation problem that 
acts as the organizing center in both the simple normal form example ([T]) 
and the Lorenz example. Specifically, we generalize the results of [IB] by 
studying an n-dimensional system of equations containing a subcritical Hopf 
bifurcation of a stable equilibrium. As in the Lorenz example, the gain ma- 
trix K for this system is such that the Pyragas feedback only acts in the 
directions tangent to center manifold of the uncontrolled system near the 
Hopf bifurcation point, and in this tangent plane the gain collapses to a 2 x 2 
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matrix that is proportional to a rotation matrix. Thus, we consider a family 
of gain matrices K that are parameterized by a magnitude Bq and a phase /3 
(c/. ([I]) written in terms of real variables). 

The additive Pyragas feedback results in a delay differential equation. We 
use methods of bifurcation theory to show that Pyragas control can stabi- 
lize the small-amplitude UPO in a neighborhood of its bifurcation provided 
60 and (3 are chosen appropriately. Specifically, our analysis applies in a 
neighborhood of a threshold value for 60, which depends on the phase angle 
P, which must lie in a particular interval that we determine. The interval 
depends only on the cubic coefficient of the Hopf normal form for the uncon- 
trolled problem. In particular, we find that this interval for /3 always exists 
provided that the imaginary part of the cubic coefficient of this normal form 
is nonzero. The threshold value for bo is associated with a highly degener- 
ate Hopf bifurcation of the zero solution of the delay differential equation. 
Specifically, for this gain modulus, the critical eigenvalue of the linearized 
problem does not cross the imaginary axis as the bifurcation parameter is 
varied so the "nonzero-speed" eigenvalue crossing condition for a generic Hopf 
bifurcation is violated. Moreover, a center manifold reduction of the delay 
differential equation to Hopf normal form reveals that the cubic coefficient 
is purely imaginary for this threshold value of 60, necessitating that one go 
to higher order than cubic in any analysis of the bifurcating periodic orbits. 
The analysis of this degenerate bifurcation problem provides the basis for our 
claims that Pyragas control can stabilize UPOs that are born from a generic 
subcritical bifurcation of a stable equilibrium in the uncontrolled problem. 
It also explains why the same sequence of bifurcations occur in the normal 
form example ([1]) as in the Lorenz example analyzed in |16] . 

The remainder of this paper is organized as follows. Section 2 reviews the 
stabilization mechanisms identified by Fiedler et al. ^ for the Hopf normal 
form example and formulates our generalized problem. Section 3 contains 
our key results. It determines the restrictions on /3 for effective stabilization. 
It also identifies and analyzes the degenerate bifurcation that acts as an 
organizing center for the control problem. Section 4 presents a rigorous center 
manifold reduction for the delay differential equation, with certain details 
relegated to an Appendix. It thereby substantiates our heuristic arguments 
made in Section 3. Section 5 summarizes our findings and discusses some 
open questions and future directions of research. 
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2 Problem Formulation 



In this section we review the mechanism of stabihzation identified by Fiedler 
et al. [9] for the Pyragas-controlled Hopf normal form. We then formulate 
the generalized problem that will be studied in this paper: an n-dimensional 
system of differential equations containing a generic subcritical Hopf bifur- 
cation of a stable equilibrium, with Pyragas-type delay terms added only in 
particular directions. We use a center manifold reduction for the uncontrolled 
problem to estimate the period of the UFO, which we use as the delay time 
r for the feedback. 

2.1 Background of Stabilization Mechanism 

Fiedler et al. [9] consider equation ([1]), where z G C and parameters A, 7 G M. 
The feedback gain is a complex number hoe^^ . For 60 = and z = re*^ we 
have 

r = (A + r^) r, 
e = 1 + 7r^ 

Bifurcating unstable periodic orbits, or "Pyragas orbits", with amplitude 
= —A, coexist with the stable trivial equilibrium for A < 0. The goal 
is to stabilize this branch of periodic orbits in a neighborhood of A = by 
adding the feedback term (69 7^ 0). The Pyragas orbits have minimal period 
r = 27r/(l — 7A), which is chosen as the delay time in ([1]). 

We now summarize the bifurcation structure associated with the z = 
solution of ([1]) in the (A, bo)-pla.ne to inform our discussion in subsequent 
sections. These results, and more details, can be found in [HI [101 US]- Fig- 
ure [m^a) shows two curves of Hopf bifurcations in the (A, 6o)-plane. One is 
the original Hopf bifurcation to the Pyragas orbit, which occurs at A = 
for every value of bo- The other Hopf bifurcation is a consequence of the 
additive delay terms and occurs along the curve 60 = ^o^(A)- It produces a 
branch of delay-induced periodic orbits, i.e. a periodic orbit that arises due to 
the addition of delay terms and one for which the feedback does not vanish. 
The two bifurcation curves intersect at the point (A, 60) = (0,&o). A curve 
of transcritical bifurcations also emanates from this point, which acts as the 
organizing center for the bifurcation structure of the problem. 

Stabilization of the Pyragas orbits involves two bifurcations. Without 
feedback, the trivial equilbrium is stable for A < and unstable for A > 0, 
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x = o 

(a) (b) 

Figure 1: (a) Schematic from [16] showing Hopf bifurcation curves A = 
(red), bo = &o^(A) (blue), and a transcritical bifurcation curve (dashed hne) 
that divide the (A,6o)-plane into five regions. The Pyragas orbit is stable 
in both the grey shaded and striped regions, (b) Schematic showing the 
bifurcations of the Pyragas and delay-induced periodic orbits as the path C 
in (a) is traversed counterclockwise about the point (A,6o) = (0,&o)- Solid 
lines indicate stable solutions and dotted lines unstable solutions. 

and the Hopf bifurcation at A = is subcritical. However, the delay terms 
can change the stability of the trivial equilibrium. For Bq > 6q, we find 
that the stability of the trivial equilibrium switches to being unstable for 
A < and stable for A > (in a neighborhood of A = 0, 6o = ^o)- Since 
both the location of the Hopf bifurcation at A = 0, and the location of the 
Pyragas orbits (in A < 0) are independent of bo, then the Hopf bifurcation 
must change criticality from subcritical to supercritical. This in turn means 
that the Pyragas orbits must now be stable. The second bifurcation involved 
in the stabilization mechanism occurs for A < 0. The Pyragas orbit is 
unstable for small values of bo, but as the feedback magnitude is increased, 
the Pyragas orbit and the delay-induced periodic orbit exchange stability in 
a transcritical bifurcation. These mechanisms can both be seen in Fig. [H 
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2.2 Generalized Problem 



The generalized problem is formulated for an n-dimensional parameterized 
system of differential equations 

y = g(y,A), 

where g : R" x R — > R" is C'^ (A; > 5, n > 2) and A is a bifurcation 
parameter. We assume there is an equilibrium solution branch y = ?/o(A) 
in a neighborhood of A = Ah, which loses stability at A = Ah as a simple 
complex conjugate pair of eigenvalues of the linear stability matrix cross the 
imaginary axis, i.e. at a Hopf bifurcation. We assume the Hopf bifurcation is 
subcritical, that is, it gives rise to a branch of unstable UPOs which coexist 
with a stable equilibrium. 

Without loss of generality we can introduce shifted variables x and bi- 
furcation parameter A so that the equilibrium is located at x = for A in 
a neighborhood of Ah, and A is proportional to (A — Ah)- Specifically, we 
define A so that, in a neighborhood of A = 0, x = is stable for A < and 
unstable for A > 0. The bifurcating branch of UPOs exists for A < 0. In 
terms of the shifted variables, the system is 

i = f(x,A), (2) 

and the Jacobian matrix -Df (0, A) has a pair of complex conjugate eigenvalues 
fj,, fi such that 

/i(A) =p(A)+za;(A), 

where 

w(0) = uo> 0, 
p(0) = 0, 

p'(0) > 0. (3) 

We further assume that the remaining n — 2 eigenvalues of Df(0,0) have 
negative real parts. 

2.3 Delay Time for Feedback 

In order to estimate the period, T, of orbits on the branch of UPOs, we 
perform an (extended) center manifold calculation [17j to reduce ([2]) to Hopf 
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normal form 

z={p (A) +iuj{X))z+ (c (A) + id (A) ) \z\h + 0{\z\^z) 
in a neighborhood of A = 0. In polar coordinates, this is 

r = p(X)r + c(A)r^ + O(r^) 

e = u{X) + d{Xy + 0{r^). (4) 

We assume that ([2]) depends smoothly on A so that the coefficients can be 
expanded in Taylor series about A = 0: 

p{X)=p^X + OiX^) 
u{X) = uo + uiX + 0(A2) 
c(A) = Co + 0{X) 

diX) = do + OiX). (5) 



Then (HI) becomes 



r = piXr + cqt^ + C(A^r, Ar^, r^) 

e = Uo + UiX + dor^ + 0{X\Xr\r^). (6) 
Neglecting the higher order terms in (jH]) we find: 

r = PiXr + CqV^ 

9 = Uq + UiX + dor"^, (7) 

and it can be shown [TS] that the dynamics of this truncated normal form 
are qualitatively unchanged when one considers the infiuence of the higher 
order terms provided pi, cq, 7^ 0. 

Pyragas orbits exist with amplitude = —piX/co for A < 0. (From ([3]) 
we have pi > 0, so for a subcritical Hopf bifurcation we must have Cq > 0.) 
These orbits have period 

^ 27r 

T ^ , (8) 

uq + (wi - P17) A ' 

where 7 = do/co captures the dependence of the oscillation frequency on 
the amplitude of oscillations. Choosing the delay such that r = T ensures 
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that the feedback vanishes when the targeted periodic orbit is reached. Our 
estimate of T in ([H]), which is based on the cubic normal form (jl]), is good 
through 0{X). 

Next, we add Pyragas-type feedback to ([2]) which gives 

x = f(x,A) + K(x(t-r)-x(t)), (9) 

where K is the constant gain matrix. As in [16], feedback is added only 
in the directions associated with the linear center eigenspace of ([2]) at the 
Hopf bifurcation point. Specifically, after a (A-dependent) coordinate trans- 
formation and a rescaling of time by the (A-dependent) delay r, takes the 
form: 




p{X) -u{X) \ /x{t)\ //l(x(t),l/(^),w(t);A)^ 

co{X) p{X) y{t) \+t\ Mx{t),y{t),w{ty,X) 

Or„, D(A) / \w(t) \Ux(t),y(t),w(t):X)j 




•(„_2,2) D(A) J \wit)J Vf,(x(t),y(t),w(t);A)^ 

bo cos f3 —bo sin f3 „ \ / x(t — 1) — x{t) ^ 



bo Sin /3 60COS/3 \ i yit -I) -yit) 

0(n-2,2) 0(„_2,n-2)/ \w(t - 1) - w(i(;) 



(10) 



where w and id are (n — 2)-dimensional column vectors, 0(ij) is a zero 
matrix with i rows and j columns, and D(A) is an (n — 2) x (n — 2) matrix in 
Jordan normal form. The eigenvalues of D(A) have negative real part for A 
sufficiently small. The feedback depends upon two parameters: an amplitude 
bo > and a phase angle /3 G [0, 27r). 



3 Hopf Normal Form for System with Feed- 
back 

We analyze the bifurcation structure of system ( |TOl) by considering the ap- 
propriate two-dimensional Hopf normal form in a neighborhood of A = 0. We 
show that as bo increases through some critical value feg, the Hopf bifurcation 
at A = changes from subcritical to supercritical, provided the phase angle 
f3 is chosen appropriately. Hence there is a range of bo > b^ for which the 
Pyragas orbit bifurcates stably. 
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We show further that the Hopf bifurcation at the point (A, Bq) = (0, 6q) 
is degenerate for two reasons: (a) the nonzero eigenvalue crossing condition 
of the Hopf bifurcation theorem is violated, and (b) the cubic coefficient of 
the normal form is purely imaginary. First we perform a linear stability 
analysis to show that the nonzero crossing condition is violated, which leads 
to an explicit formula for 6q. Another linear consideration, specifically the 
requirement that the Pyragas branch bifurcates from a stable equilibrium, 
determines restrictions on the parameters 7 and (3. In particular, we must 
assume that 7 7^ 0, and we require that /3 lie within a specified range. 

We use results from the linear analysis, together with information on 
the Pyragas branch, to argue that the real part of the cubic coefficient of 
the Hopf normal form also vanishes at A = 0, 60 = ^o- This result is later 
substantiated in Section H] (with details in Appendix B) by a center manifold 
reduction of the delay differential equation (|T0|) . The degeneracy at cubic 
order necessitates that quintic terms in the Hopf normal form be retained. 
The bifurcation analysis for this problem is performed at the end of this 
section. 

3.1 Degeneracy Associated with the Linear Normal 
Form Coefficient 

We first perform a linear stability analysis to determine the feedback magni- 
tude at which the eigenvalue crossing condition is violated. Since the feedback 
terms in flTU]) only act in the center directions, and all the other coordinates 
are linearly decaying, we focus our linear stability analysis on the x and y 
equations. In terms of the complex variable z = x + iy we obtain the linear 
delay differential equation: 

= r(A) (p(A) + tuiX)) z{t) + r(A)6oe^^ {z{t - 1) - z{t)) , 

with (3 G [0, 27r) and 60 > 0. Solutions take the form z = e^^, where rj satisfies 
the characteristic equation xiv) = 0) with 

Xiv) ^ r(A) (p(A) + tuiX)) + r(A)6oe^^ (e"^ - l) - V- (H) 

At A = 0, this equation becomes 

2m + —boe'^ (e-" - l) - 7] = 0, 
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since r(0) = 2-n/ujQ, p(0) = 0, and (X'(O) = ujq. This has a solution rj = 
rjQ = 2TTi, independent of Bq, which is as expected since the original Hopf 
bifurcation is not affected by the feedback. 

We next evaluate the eigenvalue crossing condition, that is, we consider 
to be a function of A with //(O) = 27ii, and compute Re[ri'{0)]. Specifically, 
we expand //(A) about A = as 

r/(A) = 27rz + r/iA + C(A2), (12) 

so ^^'(O) = rji. The Taylor series for p{X) and w(A) were given in ([5]). The 
delay t(A), given by ([8]), can be expanded to obtain 

r(A) = ro + riA + 0(A2), 

where 

2n 

To =— > 

Wo 

-27r (wi - pi7) 

Tl = 2 • 

Equating terms at C(A) gives 

_ 2ttpi (1 + z7) 
~ uo + 2-KhoeP ■ 

We find that the crossing condition is violated, that is. Re [//'(O)] = when 

27r (cosp + 7smp) 

In order to have a positive (finite) value for 6q, the phase angle /3 must satisfy 
the restriction 

cos/3 + 7sin/3 < 0, I3^0,tt. (14) 

Moreover, as discussed in [H], the feedback introduces additional delay- 
induced instabilities of the z = solution of ffTOj) . In order to ensure that the 
Pyragas branch can bifurcate from a stable solution, we need to ensure that 
bo is below a /3-dependent cut-off where these additional instabilities set in. 
These considerations will determine a more stringent inequality for /3, which 
can be met provided that 7 7^ 0, where 7 = d{0)/c{0) in (jl]). 



12 



We know that as Bq (i.e. the feedback vanishes), solutions of the 
characteristic equation f lllj) at A = include rjQ = 27ri, and a countable set 
that have real parts tending to —oo [19\. We determine the value of 60 (with 
A = 0) for the onset of the first delay-induced bifurcation of z = by seeking 
solutions Tj = iu of (fTTl) for A = 0. Specifically, {i^,bo) satisfy 

cos(/3 — u) = cos (3 
bo (sin(/3-z/)-sin/3) = (z/-27r)^. 

These equations have one solution at u = 2tt for all values of bo (which cor- 
responds to the original Hopf bifurcation), and a sequence (z^, 60) = {^j,bl), 
indexed by j G Z, where 

Vj = 2(3 + 27rj, K = - 



271 \ sin f3 

For the Hopf bifurcation with z/ = 27r to be from a stable equilibrium 
when 60 = we need to ensure that, for each j, either 5g > 5g or 6q < 0. 
Since /3 G (0, 27i) with /3 7^ tt, this condition is satisfied if 



where 



< 6^ < 6°, 



2^ (cos/3 + 7sin/3) ^^^^ 

^' = -2^[-^)- ^''^ 
Hence we require that /3 G (0,27r) satisfies 

sin D 

cos/3 + 7sin/3 < - — — , P^tt- (17) 

/5 — TT 

Note that < for /3 G (0,27r), so automatically ensures that I HM 

is satisfied. 

Claim: If 7 7^ 0, then there exists a value of (3 such that (fT7|) is satisfied. 
Proof: We rewrite (IT7|) as 

\ ■ iR \ ^ -sin(/3 - tt) 

— COS(/5 — TT j — 7 Sm(/5 — tt) < 



/3-7r 
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So, equivalently, we must show that there exists some /3 G (— tt, tt) \ such 
that 

~ . ~ sin/3 , . 

cos/3 + 7sm/3 > (18) 



which we can rewrite as 



+ - /3min) > (19) 

where tan/3jnin = 7- (Note, for (3 = f3mm, one obtains the smallest possible 
value of 6q given by f lT^ ). Now, observe that (sin/3)//3 < 1 for /3 G (— 7r,7r)\0. 
Since 1/1 + 7^ > 1, then it is clear that there is an open interval of /3 for 
which fll9p is satisfied. This interval will include the point where cos(/3 — 
/3min) = I5 that is, where (3 = (3min = tan~^(7). □ 

Note that when 7 = 0, condition f|T8l) becomes cos/3 > sin/3//3, or equiv- 
alently, 

tan/3</3, for /3 G (-vr, -7r/2) U (0, 7r/2), 
tan/3 >/3, for /3 G (-7r/2, 0) U (7r/2, vr), 

which is clearly never satisfied. If /3 = ±7r/2, then (sin /3)//3 = 2/77 and ( fTSl) is 
still never satisfied. Since condition (|T8l) does not hold for any /3 G {—tt, tt) \0 
when 7 = 0, the mechanism of stabilization which we discuss in this paper 
is not valid for 7 = 0. 

Henceforth we assume that 7 7^ and that /3 is chosen to satisfy f|T7|) . 
We focus on a neighborhood of A = 0, 60 = expanding the eigenvalue rj 
in ( fT2|) in a two-variable Taylor series in order to zoom in on the behavior 
near the point (A, 60) = (0,feo)- 

r/(A, 6) = 2m + 7/2^ + + V^XS + r]^5^ + 0(A^ A'^, A(5^ (5^) (20) 

where 

5 = h,- hi. 

Note that there is no linear term in A, since we already know that at 5 = 0, 
//'(A) = 0. Also, when A = 0, // = 27ri is a solution to xiv) = ^) which 

immediately implies that 772 = ^75 = 0. We assume, generically, that Re{rj^) 7^ 
0. Note that to determine 773 we would need, among other quantities, the 
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(a) (73 < (b) (73 > 



Figure 2: Schematic depiction of the stabihty of the zero equihbrium in the 
(A, (5)-plane vahd in a neighborhood of A = 5 = from (122|) when (a) < 
and when (b) cr^ > 0. Stable (unstable) regions are shaded (unshaded) and 
marked with s (u). 



delay time r through order A^, which cannot be computed using the cubic 
truncation of the Hopf normal form ([71). 

Substituting (l20i) into (fTTl) and equating terms at 0{\6) yields 



4pi7r^ (cos /3 + 7 sin Z?)' 
~i (72 + l)sin2/3 



i?e(774) = —2 , , ,^,,..2^ < 0' (21) 



where the inequality follows from (|3]) and (fT^ . Finally, defining a{X,6) = 
Re{ri), we have 

(t(A, 5) = asA^ + a4A5 + C»(A3, X^5, XS^) (22) 

where (T4 = Re{rj4) < by (pTj) . and (T3 = Re{rj^) 7^ (generically) . From 
this we are able to deduce the arrangement of the regions of stability of the 
zero equilibrium around the point {X,5) = (0,0). The arrangement depends 
on the sign of (J3, and the two cases are shown in Fig. [2l 



3.2 Degeneracy Associated with the Cubic Normal Form 
Coefficient 

The linear stability analysis reveals that the nonzero eigenvalue crossing con- 
dition is violated at bo = Bq. We now argue that the cubic coefficient of the 
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normal form equation is purely imaginary aX S = bo — Bq = 0, which is the 
second degeneracy of the Hopf bifurcation at (A, 5) = (0,0). This follows 
directly, as we now show, from the linear calculations in Section 13.11 and the 
fact that the existence of the Pyragas orbit is unaffected by the addition of 
the feedback terms. 

According to bifurcation theory for dynamical systems, including delay 
equations, flTU]) restricted to its center manifold in a neighborhood of the 
Hopf bifurcation point at A = can be converted into normal form via a 
series of nonlinear, near-identity coordinate transformations pTl [20| [211 [22] . 
This normal form is 



where z, ri,k,q G C. Here we have retained the quintic terms in anticipation 
of the result that Re{k) = at A = 5 = 0. We will assume that Re{q) ^ 0, 
and analyze the quintic truncation of f[23l) . We demonstrate for a specific 
numerical example in Appendix A that the coefficient of the fifth-order term 
does not vanish, which we expect to be true generically. Rewriting (I23p in 
polar coordinates, truncating the terms above fifth order, and considering 
only the real part of the equation, we have 



where a = cr(A,(5) is the expansion given by f[22|) and = /c/j(A,5), qn = 
gij(A, 6) are the Taylor series expansions for the real parts of the coefficients 
of the cubic and quintic terms, respectively. 

The zeros of correspond to the limit cycle solutions on the center 
manifold of the original problem f[TOl) in a neighborhood of A = 0. (Although 
the periodic orbits are not circular in the original coordinates on the center 
manifold, the normal form transformation is the coordinate transformation 
that makes them circular.) The Pyragas orbit is, by construction, unaffected 
by the control and must exist as a zero of for A < 0. From (jlj), we have 
that the Pyragas orbit satisfies = — p(A)/c(A), for which the first order 
approximation is = —piX/co- We then factor to obtain 



z = rjz + k\z\'^z + qlzl'^z + 0{\z\^z). 



(23) 



(cr + knr^ + qur^) , 



(24) 



r = 



r 



(p + cr'^) (a + 6r^) , 



(25) 
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where a is known from the hnear problem, and p, c are known from the 
original uncontrolled Hopf bifurcation and are independent of 5. 

To find a Taylor expansion of in terms of A and 5, we substitute the 
Taylor expansions of p, c (given by (|5]) ) and the Taylor series of a (given 
by ([22])) to obtain 

ka{\ 5) = '-^ + ^ + ^ + 0(A^ X5, n 
Pi Co pi 

It follows immediately that fc_R(0, 0) = 0. We have thus shown that the real 
part of the cubic coefficient of the normal form vanishes at (A, 6) = (0, 0). 

3.3 Bifurcation Analysis 

We now analyze ([5]) in a neighborhood of (A, 6) = (0, 0) with a goal of 
determining all qualitatively distinct bifurcation diagrams associated with 
the distinguished bifurcation parameter A, and showing that there is a region 
of parameter space in which the Pyragas orbit is stable. 

Expanding all coefficients in f l25|) in Taylor series in A and 6 = — b^, we 
have, at leading order, 

r = r (piA + cor^) (aiA + + qor"^) , (26) 

where 

tti = — , 
Pi 

02 = — < 0, 
Pi 
qm 

Qo = • 

Co 

The inequality for 0:2 follows from the fact that pi > and (74 < 0. 

The Pyragas orbits always exist for A < 0. On the other hand, the delay 
terms in ffTOl) create additional Hopf bifurcations that give rise to delay- 
induced periodic orbits. From f l26|) . we see that the delay-induced periodic 
orbits exist provided 

r, —a-\ A — 
= — > 0, 
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from which we find the delay-induced Hopf bifurcation fine 

S = 6oj = ^ (27) 

in the (A, (5)-plane. If go > 0, the delay-induced periodic orbits exist in the 
region of the (A, (5)-plane above this line, while if go < they exist in the 
region below this line. 

A transcritical bifurcation of periodic orbits occurs when the Pyragas 
and delay-induced periodic orbits have the same amplitude, so that the right- 
hand side of fl26|) is a perfect square. The equation for the line of transcritical 
bifurcations is 

, = ^,„=(£l*ZLi^l££)A, (28) 

with Pi,cq > 0. Note that f l27|) and fl28|) are only correct through 0{X) and 
that this bifurcation curve must terminate at A = since the Pyragas orbits 
only exist for A < 0. 

The arrangement of the curves flTTI) and fl251) in the (A, 5)-plane gives 
two qualitatively different cases for the bifurcation structure of ( l26l) . each of 
which may be divided into three subcases when the distinguished bifurcation 
parameter A is varied with fixed 6. Recall that Co,pi > and 02 < 0; we 
assume the generic situation where ai,go 7^ 0. If goai < 0, then the sign 
of pigo — aiCo is determined, but otherwise, it can be of either sign. (If 
Piqo — aiCo = 0, then the line of transcritical bifurcations occurs at 6 = 
5tc = 0.) The possible cases are thus: 



Case (1): go > 






(a) 


ai > 


Piqo - ttico > 


(b) 


ai > 


PiQo - aico < 


(c) 


ai < 


Pigo - aico > 


Case (2): go < 






(a) 


ai > 


pigo - aico < 


(b) 


«i < 


Piqo - aiCo < 


(c) 


ai < 


pigo - aico > 



(29) 



For each subcase, we show the regions of existence of the Pyragas orbits 
and delay-induced orbits, their stability, and the stability of the trivial equi- 
librium in Fig. 131 Bifurcation diagrams showing the amplitude of the Pyragas 
and delay-induced periodic orbits as a function of the parameter 6 are shown 
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(a) Case (la) (b) Case (lb) (c) Case (Ic) 




(d) Case (2a) (e) Case (2b) (f) Case (2c) 

Figure 3: Bifurcation sets for fl26p in the (A, 5)-plane for the quahtatively 
distinct cases from ( 129|) . Curves of Hopf (transcritical) bifurcations are indi- 
cated by a sohd (dotted) black hne; 6 = 6di are the delay-induced bifurcation 
curves and 6 = 6tc are the transcritical bifurcation curves. The trivial equi- 
librium is stable in the light and dark grey shaded regions. The Pyragas 
orbit is stable in the dark grey and striped regions. The delay-induced orbit 
is stable for Sdi < S < 6tc in Case 1, and for 6 < 6tc in Case 2. 

in Fig. lU In both Cases (1) and (2) the Pyragas orbits are stabilized as soon 
as 6 is increased beyond the transcritical bifurcation. However, for Case (1), 
there is a smooth transition from the stable delay-induced orbit to the sta- 
ble Pyragas orbit, whereas for Case (2), the zero solution can coexist stably 
with the delay-induced or Pyragas orbit, and hence hysteresis is expected. In 
Fig. Owe present bifurcation diagrams showing the amplitude of the periodic 
orbits vs. the distinguished bifurcation parameter A. Six distinct cases are 
manifest, but, if 5 > 0, there is always a region in which the Pyragas orbits 
are stable in a neighborhood of A = 0. It is on the basis of this observation 
that we assert that Pyragas control can stabilize the subcritical branch of 
periodic orbits, provided 7 7^ and (3 satisfies (IT7|) . 
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I|x| 



DL 




H 



(a) Case (1) (b) Case (2) 

Figure 4: Bifurcation diagrams indicating the amplitude of the periodic orbits 
vs. the parameter 5 = 6o ~ for fixed A < in fl2^ . Sohd (dotted) fines in- 
dicate stable (unstable) solutions. Tfie Hopf bifurcation to tfie delay-induced 
orbit is indicated by an H, and tfie transcritical bifurcation is indicated by 
an open dot. Tfie delay- induced periodic orbits are labeled DL Tfie Pyragas 
orbits are unlabeled red lines, wfiicfi are fiorizontal because tfieir amplitude 
is independent of 6. 

3.4 Connection with Singularity Theory 

Tfie degenerate Hopf bifurcation ( l23l) can be related to a degenerate steady- 
state bifurcation problem witfi a Z2 symmetry by focusing on tfie r equation. 
Tfiis takes tfie form 

r = s{u; A)r, u = r^, 

witfi defining condition s(0, 0) = ensuring tfiat a bifurcation of tfie r = 
solution occurs at A = 0. Tfie degenerate bifurcation of interest is defined by 

|(0.0)^|(0.0)^0. 

Tfiis is CO dimension-two as a bifurcation problem, or a codimension-tfiree 
pfienomena (i.e. including s(0, 0) = 0). Tfiis bifurcation problem is analyzed 
using metfiods of singularity tfieory in tfie book by Golubitsky and Scfiaeffer 
(Cfiapter VI) ^23j. Tfiey prove tfiat it fias tfie normal form 

r = r (air^ + 2mAr^ + a2A^) 

provided 7^ aia2, wfiere, by suitable rescaling, it is possible to set Oi = ±1 
and 02 = ±1. Tfiey analyze its universal unfolding 

r = r (flir^ + 2mAr^ + a2A^ + 04 + 2a5r^) . 
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6 < 




(a) Case (la) (b) Case (lb) (c) Case (Ic) 




(d) Case (2a) (e) Case (2b) (f) Case (2c) 

Figure 5: Bifurcation diagrams showing the amphtude of the periodic orbits 
vs. the parameter A for (|26ll (see (l29l) for description of cases). Sohd (dotted) 
hues indicate stable (unstable) solutions. Hopf bifurcations are indicated by 
an H, and the transcritical bifurcation is indicated by an open dot. The 
delay-induced periodic orbits are labeled DI, while the Pyragas orbits are 
unlabeled red curves. 

Because we assume that the time delay r coincides exactly with the period 
of the UPO, the unfolding parameter 04 = in our problem. Specifically, we 
have 05 oc 6, and we retain a term proportional to X6 in place of 04. The 
singularity theory unfolding results are expected to apply directly if we were 
to consider deviations of r from the period of the UPO. 

4 Center Manifold Reduction of the Delay 
Differential Equation 

In this section, we use a center manifold reduction for delay differential equa- 
tions to confirm that the cubic coefficient of the normal form (|23l) is purely 
imaginary at (A, 60) = (0,&c). The theory is well-developed and is described 
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thoroughly in, for example, [211 In general, because the center manifold 
cannot be determined exactly, an approximation must be constructed, and 
this calculation can be facilitated by using a computer algebra program such 
as Maple™ [26]. 

We perform the reduction at A = so that r = r(0) = 27r/a;o. Here 
we will focus on the simple case in which fllU]) has no quadratic nonlin- 
earities when Taylor-expanded about the origin. We relegate the general 
case, in which quadratic nonlinearities are also present, to Appendix B. We 
rewrite (fTO!) as 

x(t) = (Jo - Go)x(t) + Gox(t - 1) + f (x(t)), (30) 
where x(t) = (x, w)^ and the matrices Jo and Go are given by 

-Wo o 

"(2'"-^) 

0(„_2,2) 

6o cos 13 —bo sin (3 
Gq = — I 6o sin /3 Bq cos /3 

0(n-2,2) 0(„„2,ri-2) 

As in f ITU]) . 0(jj) is an z x j zero matrix and D(0) is an {n — 2) x{n — 2) matrix 
in real Jordan normal form containing the decaying eigenvalues. For the case 
we consider here, after Taylor-expanding about the origin, the vector field of 
nonlinear terms takes the form 

f (x(t)) = f3 + . . . , (31) 

where 

fs = — f2Ax' + f2,X^y + f2,Xy^ + /272/' + • • • (32) 

and (cubic) terms containing any of the Wi, i = 1 . . .n — 2 have not been ex- 
plicitly written because they do not contribute to the subsequent calculation. 
The quantities /u . . . f27 are scalars and fd4 . . . f^r are vectors of dimension 
n-2. 

We follow [22] in performing the center manifold reduction of fl30p to Hopf 
normal form at A = 0. In order to construct an appropriate phase space for 
the solutions of the delay differential equation, we define: 

^tio) = x{t + e), -i<e<o. 
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(Recall we rescaled time by the delay r = l-ajuj^^ so that the delay time is 
fixed to be equal to 1.) We write (15(1 as a functional differential equation 

x(t) = L(xO + f (Xi) 

evolving in the Banach space B = C([— 1,0] ,R") [25], where L is a linear 
mapping defined by 

= (Jo - Go)0(O) + Go0(-1), (33) 

and f is a nonlinear functional defined by 

f(0) = f(0(O)), (34) 

where G S. Here we have only 0(0) and 0(— 1) on the right-hand sides of 
(!33|) and (IMI) because there is a single fixed delay that appears only in the 
linear terms of (130!) . Linearizing (!33|) about the trivial solution x(t) = we 
obtain 

x(t) = L(xi). (35) 

As discussed in [2l], at the Hopf bifurcation point, the characteristic equation 
associated with (!35|) has two purely imaginary roots. We are interested in the 
case where the remaining (infinite number) of roots have negative real parts. 
(The conditions to ensure this led to the restriction placed on /3 determined 
in Section [31 ) Hence, the solution space of (l35l) can be decomposed as B = 
N (B S, where is a two-dimensional "neutral" (center) eigenspace spanned 
by the solutions to (l35l) corresponding to the eigenvalues with zero real part, 
and S is the (infinite-dimensional) stable eigenspace. 

The goal of the decomposition is to determine the bases needed for de- 
riving the two-dimensional equation that governs the dynamics on the center 
manifold. When system (l30!l contains only cubic nonlinearities, the equation 
can be obtained by calculating just two quantities associated with the linear 
problem, since the tangent plane approximation to the center manifold ap- 
plies at leading order. First we need a basis ^{0) for the center eigenspace 

of the linear problem, with 6 G [—1,0]. Second, we need a basis ^{^) 
for the center eigenspace of a linear problem dual to (l35l) . with ^ G [0,1]. 
The system dual to (l35l) is defined via the bilinear form given in |23]; the 
bilinear form is also used to normalize ^'(0- To summarize, the equation on 
the center manifold takes the form of a two-dimensional ordinary differential 
equation, which to cubic order is 

u = Hu + *(0)f3($(0)u) (36) 
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where 



H 



27r 
-27r 




m =\ 0-1 , (37) 



and 



Bsm/3 - {A + B cos 13) 



with 



A 



Wo 



+ 4c<;o7r6o cos /3 + Air'^bl 



B = ^^^0^0 (39) 

+ 4wo7r6o cos /3 + 47r262 

In complex form, f l36p becomes 

27r 

i = 27riz + — {A + Be-'^) F^{z, z). (40) 

Wo 

Note that the only place the delay feedback parameters appear in (HOl) is in 
the prefactor (^A + i?e~*^). The real part of the coefficient of \z^z in F3, 
which is the only term that survives the normal form transformation [T7], 
takes the form 

A;^ = (A + S cos /3) Co + (5 sin /3) rfo (41) 



=co 

where 7 = (io/co, and 



A + 5(cos/3 + 7sin/3) 



3 113 
Co = -/i4 + -/16 + g/25 + g/27 (42) 

C^O = ~g^l7 — g/l5 + g/26 + g/24- 

The quantities co and in ( l42l) are the real and imaginary parts, respectively, 
of the cubic coefficient in the uncontrolled problem restricted to the center 
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manifold evaluated at A = 0. We can see that kji = when A + B{cosf3 + 
■J sin (3) = 0. Using ( 15^ we find that the value of bo for which this occurs is 
when 

, _ , c _ 

°~ ° ~ 27r(cos/3 + 7sin/3)' 

Using center manifold theory, we have thus shown that the cubic coefficient 
of the normal form vanishes at 6o = ^O) which confirms the results based on 
the heuristic argument given in Section 3. 



5 Discussion 

We have shown that the branch of small amplitude UPOs created in a generic 
subcritical Hopf bifurcation from a stable equilibrium can be stabilized using 
an appropriate Pyragas-type feedback, provided that 7 7^ and (3 satis- 
fies f|T7|) . where 7 is proportional to the imaginary part of the cubic coeffi- 
cient of the Hopf normal form for the uncontrolled problem. Specifically, we 
followed [TB] and considered feedback only in the directions that are tangent 
to the center manifold at the bifurcation point. As in [9J, the feedback gain 
matrix is parameterized by a gain amplitude, bo, and a phase angle, /3. This 
choice of gain matrix reduces the problem of choosing components of a 
matrix, to that of choosing just the two parameters 60 and /3. Moreover, 
we note that knowledge of 7 is already assumed for Pyragas control since it 
is needed to determine the time-delay for the feedback at leading order in 
A. No other information about the nonlinearities of the problem is required, 
making this a particularly simple and elegant control. 

The behavior of the controlled system for bo near 6q is governed by a highly 
degenerate Hopf bifurcation problem. In particular, both of the nondegen- 
eracy conditions for a generic Hopf normal form are violated simultaneously 
: (a) the eigenvalues of the linearized problem do not cross the imaginary 
the bifurcation parameter is varied, and (b) the real part of the cubic 
coefficient of the Hopf normal form, which determines whether the bifurca- 
tion is subcritical or supercritical, vanishes. We performed a center manifold 
reduction of the governing delay differential equations to prove this result in 
Section |H The reason this degeneracy may be surprising is that only two 
parameters (A and 60) are varied, and yet a codimension-three problem is 
obtained. In Section |3] we demonstrated that this is due to the fact that the 
time-delay of the Pyragas feedback is fixed to coincide with the period of the 
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original UPO. A consequence of this restriction is that the existence proper- 
ties for the bifurcating branch of targeted periodic orbits is unaffected by the 
feedback, although its stability is altered. The transcritical bifurcation from 
classic bifurcation theory provides an apt analogy. This bifurcation typically 
occurs when there is an (equilibrium) solution whose existence is unaffected 
by the value of the control parameter, e.g. some trivial equilibrium state 
X = that exists for all parameter values. For problems with such struc- 
ture, the transcritical bifurcation supplants the saddle-node bifurcation as 
the generic steady state bifurcation problem [23] . 

Just as one can unfold a transcritical bifurcation with imperfections, our 
degenerate Hopf bifurcation problem at (A, 60) = (0, &o) could be unfolded 
further by considering detunings of the delay time r from the period of the 
UPO. (This unfolding would also allow us to apply directly more of the 
singularity theory results for the degenerate bifurcation problem, developed 
by Golubitsky and Schaeffer [23] and described in Section |3l) We note that 
Just et al. [To] have investigated bifurcations with varying r in the setting 
of the Pyragas controlled subcritical Hopf normal form ([T]). They analyze a 
series of bifurcations in the three-parameter space r-A-60 (in our notation). 
Since they allow r to vary, the point A = 0, r = 27r, 60 = &o in their case is a 
codimension-three point, and there is a complicated sequence of bifurcations 
about this point. 

The fundamental mechanism at work for the Pyragas control of UPOs 
arising from a subcritical Hopf bifurcation was first explained in the paper of 
Fiedler, et al. [9]. They showed that the feedback introduces delay-induced 
instabilities of the bifurcating equilibrium, thereby altering its stability prop- 
erties in a neighborhood of the original Hopf bifurcation. In this way, the 
equilibrium can change from being stable for A < to being stable for A > 0, 
and then the original subcritical Hopf bifurcation (with increasing A) can be 
converted to a supercritical Hopf bifurcation (with decreasing A). This under- 
standing suggests a number of interesting directions for generalizing Pyragas 
feedback for other bifurcations, where a similar simple control mechanism 
may work. 

The Hopf bifurcation is just one of several bifurcation mechanisms by 
which UPOs can generically arise in dynamical systems. As mentioned in the 
introduction, two other generic mechanisms for creating UPOs in dynamical 
systems include homoclinic (or heteroclinic) global bifurcations and saddle- 
node bifurcations of limit cycles. The application of Pyragas-type feedback 
to stabilize UPOs arising in some of these situations has already been studied 
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in the context of specific examples [3 |H], and in both the saddle- node and 
heteroclinic case, the UPO is stabilized in a steady-state bifurcation with a 
delay-induced periodic orbit. In the heteroclinic case, the original problem 
was again of a higher dimension than the UPO. The gain matrix was chosen 
as in this paper, to be tangent to the two-dimensional manifold containing 
the UPO. We emphasise that this approach, studying the generic bifurcations 
which generate UPOs, leads to results which allow one to choose a gain matrix 
and predict the parameter region for which control can be achieved, with 
very little knowledge of the structure of the UPO. Preliminary calculations 
attempting to stabilize a UPO arising from a generic homoclinic bifurcation 
have indicated that this is in fact not possible, in distinct contrast to the 
other cases P?] . 

An additional promising direction for generalizing Pyragas feedback for 
control of UPOs is to Hopf bifurcation problems with symmetry. A feature 
of such bifurcations, which arise naturally in pattern forming systems [28], is 
that they typically involve multiple complex conjugate pairs of eigenvalues 
crossing the imaginary axis simultaneously. Thus the results of this paper 
do not apply directly since we assumed that a simple complex conjugate pair 
crosses the imaginary axis at A = 0. Moreover, in equivariant Hopf bifurca- 
tion problems, a number of periodic orbits, distinguishable by their group of 
symmetry properties, may bifurcate simultaneously [29]. This group theo- 
retic classification result suggests investigating feedback controls, analogous 
to Pyragas, that exploit the targeted symmetry properties (which may be 
spatial, temporal, or a combination of spatio-temporal symmetries) in such 
a way that the feedback vanishes on the targeted state only. The goal would 
be to use the feedback control to change stability properties of the bifur- 
cating equilbrium, while naturally preserving the bifurcation to the targeted 
state. Some first steps in this direction of controlling patterns arising in 
a Hopf bifurcation, using extensions of Pyragas control tailored to spatial 
symmetries (as well as temporal ones), were shown to work for certain plane 
wave solutions of the complex Ginzburg-Landau equation in regions of the 
Benjamin- Feir unstable regime [5D| 1^. Further work on the stabilization of 
unstable patterns using these methods is ongoing. 
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Appendix A 

Our analysis of the normal form for the degenerate Hopf bifurcation prob- 
lem fl2^ associated with A = 0, 60 = ^q, presented in Section 3, is valid 
provided that certain nondegeneracy conditions are met. For instance we re- 
quire that (T3 7^ in (!22l) so that the coefficient of the linear term has real part 
proportional to at 60 = ^o- Moreover, we assume that at (A, 60) = (0,&o), 
the real part of the coefficient of the quintic term in the normal form is 
nonzero. We do not perform a center manifold reduction of the delay differ- 
ential equation to fifth order to determine this quintic coefficient since our 
results show that the important qualitative features are the same whether 
the quintic coefficient is positive or negative. In order to allay any possi- 
ble concerns that there is a hidden structure to the problem associated with 
Pyragas control that would force the real part of the quintic coefficient to 
also be zero at (A, 60) = (0,&o), we present a simple numerical example here 
for which we show that the coefficient is nonzero. We expect therefore that 
this nondegeneracy condition will be satisfied generically. 

Our numerical example is based on the Hopf normal form equation ([1]) 
with Pyaragas feedback control that was investigated by Fiedler et al. [S]. 
After rescaling time by r so that the delay is fixed, we have 

z{t) = r(A + i)z{t) + r(l + i^)\z{t)\^ z{t) + Th^e'\z{t - 1) - z{t)) (A.l) 

with r = 27r/(l — 7A). For our numerical example, we choose parameter 
values 7 = —10, (3 = tt/8. Following the procedure in Sectional a linear 



28 



-4 



X 

X 



-10 -9 -8 -7 -6 -5 -4 -3 -2 

Figure 1: The figure shows log(r) vs. log |A — Xdi\, indicated by x's, for the 
delay-induced periodic orbit computed from flA.l|) using dde-biftool [32] 
for 6 ~ —0.0012. For comparison, a line with slope 0.5 is plotted above the 
data points. 

stability analysis yields 

6^ ^ 0.055 
^3 ^ -247 
(74 ~ —65, 

where the aj appear in the expansion of the linear growth rate a = (TsA^ + 
(J4A5 + ■ ■ ■ associated with perturbations of the zero solution. 

Since as 7^ 0, the normal form can, from ( l26ll . be written in the following 
factored form, valid in a neighborhood of (A, Bq) = (0, 6q): 

r = r(A + r^)(o-3A + (746 + qmr^). 

We want to show that qro 7^ 0, so we focus on the delay-induced periodic 
orbits with amplitude r satisfying 

asA + (74(5 + QRor^ = 0. (A.2) 

For a fixed 6, the value of A at which a Hopf bifurcation to a delay-induced 
periodic orbit occurs, which we denote by Xdi{S), is defined by 

o-3Ad/(5) + (^aS = 0. 
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Then (1A.2P becomes 



(T-s (A - \di{S)) + qRor^ = 0, 

where 

X^j{S) = ^ ^ -0.264 S. 

-as 

We followed the branch of delay-induced periodic orbits in a neighbor- 
hood of \ = Xdi for a sequence of values of 6 approaching 6 = using the 
bifurcation package dde-biftool [32]. If g^jo 7^ and A sufficiently close to 
Xnj, then we expect the amplitude r of the delay-induced branch to scale as 
\/J)^^-^^^jji\. In Figured] we graph log(r) vs. log |A — Xnil for 6 ~ —0.0012, 
and find the square root dependence of r on |A — A^/l, which is consistent 
with g^jo 7^ 0. If qro = 0, then we would not expect this scahng. 



Appendix B 

This Appendix presents the center manifold reduction for the case in which 
system flUI]) contains quadratic nonlinearities, as well as cubic. The goal is the 
same: to determine the two-dimensional equation that governs the dynamics 
on the center manifold. However, because the tangent plane approximation 
to the center manifold no longer applies, we must now approximate the center 
manifold. 

We begin by rewriting (fTOl) in the form ( l30l) . The vector field of nonlinear 
terms can be expanded as: 

f (X(t)) = f2 + f3 + . . . 

where fs is given by (1321) and 

2^ / fn^^ + fi2xy + /i3l/^ + Er=i^ (giiXWi + kuywi) \ 

f 2 = — /2lX^ + f22Xy + /232/^ + YTi=l {92iXWi + k2iyWi) \ + ... 

\ fdix^ + fd2xy + id-sy"^ + YJiZi (gdiXWi + kdiywi) J 

Here we have neglected to write terms quadratic in Wi, since they don't enter 
the subsequent calculation. 

With the bases ^{9) and ^'(6') already given by fl37j) and fl38l) . respec- 
tively, the final quantity to be calculated is an approximation to the center 
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manifold, which will be denoted by h(u, 6). To obtain terms up to cubic or- 
der for the evolution equation on the center manifold, we need only consider 
terms up to quadratic order in the equation for h(u, 6), that is, we write 

/ hi_ui6)ul + hi_i2{e)uiU2 + hi_22{e)ul 

h2(^, U) = h2_ll{e)ul + h2_l2{9)uiU2 + h2_22{0)ul 
\ hd_u{6)ul + hd_l2{0)UiU2 + hd_22{d)ul 

where hi_ii{9), hi_i2{0), . . . h2_22{0) are scalar functions of G [0, 1] and hd_ii{9), 
hd_i2{0), and hd_22{0) are (n — 2)-dimensional vector- valued functions. As 
shown in [26], h2 must satisfy the equation 



de du 
subject to the boundary condition 

dh2 



u)Hu + $(^)*(O)f2($(0)u) 



du 



Hu + $(O)*(O)f2($(0)u) = L(h2(0, u)) + f2($(^)u) 



(B.l) 



(B.2) 



where H is given by ( 137|) . Solving f IB.ip . ( ]B.2p for the coefficients of the first 
two rows of h2 yields 



h2(0) 



3^0 



1- A- Bcosl3)l3 -(Ssin/3)l3 
(5sin/3)l3 {1 - A - B cos (3)h 



(B.3) 



where I3 is the 3x3 identity matrix. A, B are given by ( 139|) and 

/ — /12 — /21 — 2/23 \ 
-2/n + 2/13 - /22 

/12 — 2/21 — /23 
/ll + 2/13 — /22 
/l2 — 2/21 + 2/23 
V 2/n + /i3 + /22 / 



(B.4) 



To determine the vectors hrf_ii(0), h^_i2(0), and hrf_22(0) we solve the following 
system of equations: 



M 



hd.ii(0) 
hd.i2(0) 

hd.22(0) 
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where 

/ -D(0) -CUoln-2 On 2 

M= 2uJoln-2 -D(0) -2cJoIn-2 

\ On 2 UJoIn~2 "0(0) 

is a matrix of size 3(n — 2) x 3(n — 2), D(0) is the real (n — 2) x (n — 2) 
matrix defined in ( !T0|) and evaluated at A = 0, 0„_2 is a (ri — 2) x (n — 2) 
zero matrix, and I„_2 is the (ri — 2) x (n — 2) identity matrix. 
Finally, the equation on the center manifold is given by 

u = Hu + *(0) (f2 ($(0)u) + Vf2 ($(0)u) h2(0) + fa ($(0)u) ) , 

or in complex form, 

271- 

z = 2niz + — [A + Be-'^) ^2.3(2, z), 

UJo 

where -^2,3 (-2) contains both quadratic and cubic terms proportional to z"^, 



etc. After performing a near identity transformation [T7], the 



real part of the cubic coefficient \z\'^z takes the form 

kR = {A + B cos (3) Co + {B sin (3) do 

= Co 



(B.5) 



with 7 
when 



A + 5(cos/3 + 7sin/3) , 
do/cQ. Once again, from equation (IB. 5 II we can see that kji 



h = hl 



27r(cos /3 + 7 sin /3) 



For completeness, we find 

3 



Co — |/l4 + |/l6 + I/25 + I/27 + 4^/13/23 — 4^/11/21 
+ 8^/1^/12 + 8^/12/13 - 8^/21/22 - 8^/23/22 

+1 £r=i^fi'H^di.ii(0) - I X]r=i^ (^^^02(0) - A;2i/idoi(0)) 
+ 8 2^i=i '^fi^i* 



^2i^a 

hdij22{^) - hihdi.uiO)) + I T.i=i k2ihdi.22i0) 



I/24 + I/26 



15 



1/17 + 



1 



1 



■(/: 



12 



1 



■(/: 



24(^0^'^^^^ ' 24aJo 

7/21/23 — T^h^flS — 



+ 



1 



■/13/: 



22 



24(^0 •^^^•^^^ 24aJo" 
i(/22)^ + ^J\if 22 



■6^(/23)^ - 6i;(/ll)^ + 2£;^/ll/22 + 2^^/23/12 
\n— 2 7 I 1 v^n-2 



+1 Ya=i 92ihdi.n{0) + I Y17=i i9iihdiA2{0) - kiihdi_n{0)) 

+ 1 Z]r=l^ (fl'2i^di.22(0) - A;2i/idO2(0)) - I X]r=l^ k2ihdi.22i0), 
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which are the precisely the coefficients c{0),d{0) in the uncontrolled Hopf 
normal form (jl]). 
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